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ABSTRACT 



Aims. 3C 345 was recently identified as a y-ray emitter, based on the first 20 months of FermZ-LAT data and optical monitoring. 
In this paper, a connection between the y-ray and optical variability of 3C 345 and properties of its parsec-scale radio emission is 
investigated. 

Methods. The Fermi-hAT data of 3C345, covering an energy range of 0. 1-300 GeV, were combined with 32 Very Long Baseline 
Array observations of the object made at 43.2 GHz in the period of January 2008 - March 2010. 

Results. The VLBA data reveal morphology and kinematics of the flow on scales of up to a 5 milliarcseconds (deprojected linear 
distances of 380 parsecs). The brightness temperature, T],(r), measured along the jet first decreases with distance oc ^-(0-95±o.69) ^^^^ 
later exhibits a break at « 0.3 milliarcseconds (mas), with Tb(r) oc ^-(4 ii ±o.85) j^j-gg^ separations. Variations of the y-ray, optical 
and parsec-scale radio emission show a similar long-term trend persistent during the entire VLBA monitoring period. The y-ray and 
optical variations on shorter time scales are related to structural changes in the jet on scales of a 0.3 mas (a 23 parsecs, deprojected), 
with the y-ray and optical flares possibly related to the evolution of four distinct superluminal components identified in the flow. 
Conclusions. The observations indicate that both the quiescent and flaring components of the y-ray emission are produced in a region 
of the jet of ~ 23 pc in extent. This region may mark the Compton-loss dominated zone of the flow and its large extent may favor the 
synchrotron self-Compton mechanism for y-ray production in the relativistic jet of the quasar 3C 345. 
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1. Introduction 

Observations of high-energy emission (keV - TeV range) pro- 
vide important insights into the physical properties of Active 
Galactic Nuclei ( AGN), in particular for the subgroup of quasars. 
However, high energy production scenarios ar e still heavily 
debat e d with many possib i lities discussed (e. g., Bottc her et al.l 
l2009t iDermer et al] 120091: iBlazeiowski et all 12004) . Most of 
these physical models depend on the emission site, mainly its 
distance to the central engine. A connection between cm-/mm- 
ra dio and high energy emission o f quasars was already suggested 
bv lValtaoia & Teraesrantal ( ll996h . |jorstad et al.l ('2001 ) and refer- 
ences therein. To investigate this connection, variability of AGN 
are studied in coordinated multi-wavelength observation cam- 
paigns covering radio to y-ray wavelengths. Their results, espe- 
cially on archetypal sources, are essential for the physical un- 
derstanding o f the emission p rocesses in AGN up to the highest 
energies (cf. lLongaiill20Tll) . These observations, in particular 
simultaneous high-resolution very long baseline interferometry 
(VLBI) observations at mm-wavelengths, provide new limits on 
the parameter space of proposed emission models, placing con- 
straints on size, energetics and locations of emission regions. 

The archetypal quasar 3C 345 is one of the best stud- 
ied "superluminal" radio sources, with its parsec-scale radio 



emission monitor ed over the past 30 years, in particular b 
VLBI (e.^ . : Birett a et al.ll98dlBaath et al.ll992l:lLobanovll99L. 
Zensus et al., ,19951: iRos eTail l2000t iKlard l2004t iJorstad et al 
2007h . Based on previous observations, spe ctral index and 



turnover frequency distributions were obtained dLobanovll 19961 
ii998a: Ros et al,, 2Q0Q,), the spec tral evolution of the jet was 
studied dLobanov & ZensusI 19991). properti e s of radio and X-ray 
emission were related (Unwinetal ] 119971: iLobanov & Rolandl 
2003l). the dynamics of the central region was investigated 
( ILobanov & Rolandl l2005h and opacity in the nuclear region 
was used to determine the ph ysical properties and matter coni - 
position of the compact jet dLobanov! 1199 8bl : lHirotanill2005h . 
Quasi-peri odic var iations of emission have been detected in 
the op tical dBabadz hanvants & Belo konlll984tlki dger & Takal 
19901) and radio dAlleretal.1 ll996riTeraesranta et al.1 1 199 



Lobanov & Zensuslll999h " regimes with a possible periodicity of 
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3.5-4.5 years and major flares occurring every 8-10 years. A 
new cycle of such enhanced nuclear activity began in early 2008 
(Larionov et al. 200g). 

The identification of 3C 345 as a y -ray source was unclear 
during the EGRET-era (1991-2000; iThompson et"anil993h 
with a possible weak detecti on between April and May, 1996 
dCasandiian &Grenieill2008l) . The y-ray emission detected by 
the Large Area Telescope (LAT), on board the Fermi satel- 
lite, in the vicinity of 3C 345 was initially associated with an- 
other quasar in the region, B3 1640-1-396 with low confidence, 
based on the 3-month data (August - October, 2008) collected 
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(lAbdo et al.ll2009h . However, with 20 months (August 2008 - 
April 2010) of LAT monitoring data, 3C 345 was identified as a 
y-ray source at GeV energies, based on multi-wa velength coun- 
terpar t locaUzations and correlated variability dSchinzel et alJ 
l201lh . 

This paper presents results from a coordinated observational 
campaign targeting 3C 345 and comb ining the 20 -month Fermi- 
LAT y-ray monitoring data (Schinz el et al.ll20ril) with monthly 
VLBI observations made at 43.2 GHz (7 mm wavelength) at the 
Section 2 discusses the observations and data reduction 
methods applied to obtain calibrated datasets for further analy- 
sis. It discusses in particular the criteria that were developed in 
order to determine the statistical significance of the 2D-Gaussian 
models (modelfits) applied to the observed radio brightness dis- 
tribution. In Section 3 we present kinematics of the pc-scale ra- 
dio jet and flux density evolution, and a possible connection with 
the high energy y-ray emission. Section 4 discusses the findings 
of this paper and in Section 5 conclusions are made from the 
results presented here. 

Throughout this paper a flat ACDM cosmology is assumed, 
with //() - 71kms~' Mpc~' and Qm - 0.27. At the redshift 
z = 0.593 (iMarziani et al.1 Il996h of 3C 345 this relates to a 
luminosity distance Z)l - 3.47 Gpc, a linear scale of 6.64 pc 
mas"' and a proper motion scale of I mas year ' corresponding 
to 34.5 c. 




-1 -2 -3 -4 

Relative Right Ascension (mas) 

Fig. 1. VLBA image of the total brightness distribution of 
3C 345 at 43.2 GHz made from observations on Nov. 30, 2009. 
Open crossed circles show the FWHM of eight Gaussian com- 
ponents applied to fit the structure observed. The shaded ellipse, 
in the lower left corner, represents the FWHM of the restor- 
ing beam. The image peak flux density is 2.1Jy beam"' and the 
RMS noise is 1 mJy beam"'. The contour levels are (-0.15, 0.15, 
0.3,0.6, 1.2,2.4, 4.8,9.6, 19, 38, 77) % of the peak flux density. 
Here the nuclear region is modeled by two circular Gaussian 
components (QO, Qll) identified as the best representation of 
the observed brightness distribution. 



2. Observations & Data Analysis 

2.1. Fermi-LAT 

The Fermi-LAJ' ( lAtwood et al.ll20()9h is a pair conversion tele- 
scope designed to cover the energy band from 20 Me V to greater 
than 300 GeV. It is the product of an international collaboration 
between NASA and DOE in the U.S. and many scientific insti- 
tutions across France, Italy, Japan, and Sweden. 

The y-ray emission of 3C 345 was identified based on cor- 
relations found between the optical variability and major y-ray 
events observed by Fermi LAT between August 2008 and April 
2010. The y-ray counterpart of 3C345 was localized to R.A. 
16''43"0.24\ Dec. -H39°48^22.7" (ISchinzel et al.ll201 ih . For this 
paper a light curve for which the y-ray monitoring data was split 
into regular time intervals, each integrating over periods of 7 
days and an energy range of 0.1 -300 GeV , was obtained in the 
fashion described in ISchinzel et al.l (1201 ih . The position of the 
y-ray counterpart was fixed to the radio localization of 3C 345. 
For the spectral shape of the y-ray emission of 3C 345 a power- 
law was used with the spectral index fixed to its 20 month av- 
erage value of F = 2.450. The particular time binning of 7 days 
in this case provides the best trade-off between time resolution 
and signal to noise. This yielded a light curve with 8 1 significant 
detections and five 2cr upper limit time intervals (JD 2454756, 
2454826, 2454903, 2454945, 2455155), in total this covers a 
time period of 602 days (20 months). In order to homogenize the 
light curve, 2cr upper limits were used as values with their error 
estimate for that interval replaced with half the difference be- 
tween that upper limit and its value determined through the un- 
binned spectral likelihood analysis. This method was app l ied for 
the calculation of the variability index in ISchinzel et al.l (1201 ih 
and lAbdo et all (1201 Ol) . 



' Very Long Baseline Array of the National Radio Astronomy 
Observatory, Socorro, USA 

^ The photon spectral index F is defined as N(E)ocE"'^, where N(E) 
is the y-ray photon flux as a function of energy E. 



2.2. Very Long Baseline Array (VLBA) 

Following the onset of a new period of flaring activity in 2008, a 
dedicated monthly monitoring campaign was initiated, using the 
VLBA to monitor the radio emission of 3C345 at 43.2, 23.8, 
and 15.4 GHz (VLBA project codes: BS193, BS194). In this 
paper only the 43.2 GHz observations are discussed, while the 
analysis of 15.4 and 23.8 GHz data is continued. The observa- 
tions were made with a bandwidth of 32 MHz (total recording bit 
rate 256 Mbits s"'). A total of 12 VLBA observations were com- 
pleted, with about 4.5 hours at 43.2 GHz spent on 3C 345 during 
each observation. Scans on 3C345 were interleaved with ob- 
servations of J13104-3233 (amplitude check, EVPA calibrator), 
J1407-H2827 (D-term calibrator), and 3C279 (amplitude check, 
EVPA calibrator). The VLBA data were correlated at the NRAO 
VLBA hardware correlator and starting from December 2009 the 
new VLBA-DiFX correlator was employed. Analysis was done 
with NRAO's Astron omical Image Process ing System (AIPS) 
and Caltech's Difmap (IShepherd et al.|[T995h software for imag- 
ing and modeling. Corrections were applied for the parallactic 
angle and for Earth's orientation parameters used by the VLBA 
correlator Fringe fitting was used to calibrate the observations 
for group delay and phase rate. A summary of all the observa- 
tions is presented in Table [T] Here, the data from the 12 epochs 
of this dedicated monitoring campaign are presented, combined 
with 20 VLBA observations from the blazar monitoring program 
of Marscher et al. (VLBA project codes BM256, SI 136) avail- 
able onlineQ. The combined data (see Table [T]i cover a period 
from January 2008 to March 2010, with observations spaced 
roughly at monthly intervals or shorter 

The brightness distribution of the radio emission was mod- 
elfitted by multiple Gaussian components providing positions, 
flux densities and sizes of distinct emitting regions in the 
jet. Fig. [T] illustrates the observed radio structure and the 
Gaussian modelfit representation. In the following, we interpret 
the eastern-most Gaussian modelfit component obtained from 

^ http://www.bu.edu/blazars/VLBAproject.html 



Schinzel et al.: Relativistic Outflow Drives y-Ray Emission in 3C 345 



3 



Table 1. Summary of VLBA Observations. 



Date 


S tot 


D 


Beam (bpa) 


Ref. 




[Jy] 




[mas] X [mas] (°) 




2008-01-17 


1.90 


2300 


0.31x0.19 (-27.2) 


1 


2008-02-29 


2.02 


1600 


0.37x0.21 (-29.1) 


1 


2008-06-12 


2.44 


2100 


0.38x0.16 (-28.2) 


1 


2008-07-06 


2.23 


800 


0.33x0.15 (-16.6) 


1 


2008-08-16 


3.78 


2700 


0.41x0.19 (-30.2) 


1 


2008-09-10 


3.67 


4300 


0.37x0.19 (-32.5) 


1 


2008-11-16 


4.43 


300 


0.37x0.33 (-1.49) 


1 


2008-12-21+ 


2.91 


5500 


0.39x0.17 (-16.9) 


1 


2009-01-24 


5.04 


8900 


0.31x0.17 (-20.5) 


1 


2009-02-19"! 


3.55 


2200 


0.37x0.20 (-19.0) 


2 


2009-02-22 


4.63 


8800 


0.35x0.15 (-19.1) 


1 


2009-03-16 


6.01 


2400 


0.43x0.30 (9.31) 


2 


2009-04-01 


5.67 


8000 


0.33x0.16 (-18.7) 


1 


2009-04-21 


5.78 


2600 


0.33x0.22 (-16.7) 


2 


2009-05-27 


7.00 


3500 


0.33x0.18 (-15.2) 


2 


2009-05-30 


8.65 


6900 


0.32x0.16 (-19.8) 


1 


2009-06-21 


5.04 


6500 


0.28x0.16 (-10.9) 


1 


2009-06-29 


6.43 


2600 


0.29x0.16 (-11.6) 


2 


2009-07-27 


7.63 


2900 


0.38x0.16 (-30.8) 


2 


2009-08-16 


5.80 


4300 


0.30x0.16 (-20.4) 


1 


2009-08-26 


6.58 


2000 


0.32x0.17 (-24.8) 


2 


2009-09-16 


6.75 


4600 


0.35x0.19(18.3) 


1 


2009-10-01 


6.67 


2300 


0.22x0.18 (-21.6) 


2 


2009-10-16 


6.80 


4800 


0.41x0.17 (-30.5) 


1 


2009-11-07 


6.38 


1400 


0.32x0.16 (-15.3) 


2 


2009-11-28 


5.31 


4900 


0.28x0.16 (-24.6) 


1 


2009-11-30 


5.53 


2000 


0.29x0.16 (-16.9) 


2 


2009-12-28 


4.95 


2400 


0.31x0.16 (-8.14) 


2 


2010-01-10 


4.25 


2400 


0.33x0.20 (-8.75) 


1 


2010-02-11 


5.57 


2100 


0.27x0.15 (-17.0) 


1 


2010-02-15 


5.11 


1300 


0.31x0.20 (-6.37) 


2 


2010-03-06 


5.38 


3300 


0.39x0.28 (-12.1) 


1 



Notes: 5,ot - total flux density recovered in VLBA image; 

D - dynamic range measured as a ratio of the image peak flux density 

to the RMS noise; 

Beam (bpa) - beam size, major axis vs minor axis with position angle 
of ellipse in parentheses; 

References: 1 ~ blazar monitoring Marscher et al. (VLBA project 
codes BM256, S1136); 2 - dedicated monitoring (VLBA project codes 
BS193, BS194). 

t - not used for the flux density analysis due to gain calibration 
problems and bad weather conditions on some of the VLBA antennas. 



the VLBI map, hereafter labeled as QO (see Fig. [U, as the base 
(or "core") of the radio jet at 43.2 GHz. The remaining fea- 
tures can signify perturbations or shocks developing in the jet. 
Locations and proper motions of these jet features are then de- 
termined with respect to QO. 

To find the best description of the observed brightness dis- 
tribution, compact emission in the nuclear region (< 0.15 mas 
from QO) was modeled using four different approaches: 1) single 
circular Gaussian, 2) single elliptical Gaussian, 3) two circular 
Gaussians, 4) three circular Gaussians. In the following discus- 
sion, these modelfitting approaches are designated IC, IE, 2C 
and 3C, respectively. 

The criterion for significant improvement between two mod- 
elfits was developed using a ratio of Teducedx^ parameters 



where the subscripts refer to the two modelfits compared. The 
order is defined by the number of degrees of freedom n of the 
individual modelfits with the condition «i < «2- The parame- 
ter xodii is obtained by solving the equation of the cumulative 
distribution function (CDF) of the distribution. 



0.68 = 



1 



r(!) 



'(2'—)' 



(2) 



where F denotes Euler's Gamma function and y represents the 
lower incomplete Gamma function. The number of degrees of 
freedom, n, depends on the number of modelfit parameters, since 
= X^^^ + Xa (n, a), with a describing the significance level 
(a - 0.68, in our ca se). The pro of and application of this prop- 
erty are discussed in lAvnild 19761) . 

Fig.|2]compares the;^'^ parameters of each modelfit, normal- 
ized to the;^'^ values of the IC (top panel), IE (middle panel) and 
2C (bottom panel). The horizontal lines in the plots represent the 
upper thresholds for significant improvement over the labeled 
Gaussian model. The IC models compared to the IE and 2C 
models have the respective threshold values of 0.77 and 0.82. IE 
models compared to 2C models have a threshold value of 0.91, 
and 2C compared to 3C models needs to have an improvement 
of 0.86 or less. For some of the observations, the representation 
of the nuclear region by IE is comparable to 2C. However, in 
order to form a consistent dataset, we adopt the representation 
by 2C providing the optimal ratio of x^ to the number of model 
parameters for all epochs. 

The errors of the obtained modelfit parameters were esti- 
mated from the im age plane, follow ing analytical appro xima- 
tions introduced bv lFomalon^ (Il999l) and iLobanovl (l2005h with 
modifications to account for the strong side-lobe case inherent 
to VLBI observations. For each component in the model, the 
root-mean-square (RMS) noise (o-p) was determined measur- 
ing around its position in the residual image. Then the respec- 
tive component was removed from the Gaussian model and the 
modified model was subtracted from the data, yielding an image 
that contained only the contribution from the component investi- 
gated. The flux density (Sp) at the peak position of the Gaussian 
was taken. Further input values were, observation beam size 
(b 



maj 



' ^min), component total flux density (St), compo- 
nent distance (r), position angle ((D), and size (d). The errors for 
the model fit parameters were calculated in the following way: 



SNR = — (SNR > 1) 



a-i 




= I ■ f ln(2.0)...ln(J^; 



(3) 



(4) 



(5) 



a-d 



p ' "lim 
5n 



= 0.5 ■ (Td 



if (c/iini > d) else (dnm = d) (6) 

(7) 
(8) 



X\ Jf0.68;l 



/0"r 180 

(1) if (r > 0) : (To = arctan — ■ 

\ r n 



(9) 
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Fig. 2. Reduced parameters of the modelfit representations of 
the core region. The top panel shows the resuhing values 
of IE, 2C and 3C Gaussian models normalized to the value of 
the IC Gaussian representation. The two horizontal lines mark 
the limit below which, for the corresponding model, a significant 
improvement (68% confidence) over IC is given (0.77, 0.82). 
Similarly the middle panel shows the X^educed values for 2C and 
3C normalized to that of IE. The horizontal line marks the sig- 
nificant improvement level compared to 2C (0.91). The bottom 
panel normalizes the x~ values to that of 2C. The dashed line 
plots the significant improvement limit for 3C (0.86). 



These error estimates only reflect the statistical image errors and 
are assumed uncorrelated. No additional systematic errors were 
taken into account, such as errors of the amplitude calibration 
that could add an additional uncertainty to the determined flux 
density values. 

3. Results 

3.1. Evolution of the radio emission in the nuclear region 

Comparisons of the component positions, flux densities and 
sizes determined from the Gaussian modelfits revealed a new 
moving emission region, labeled Q9, first detected in the VLBA 
image from June 16, 2008, followed by detections of another 
new component on January 24, 2009 (QIO), a third one on July 
27, 2009 (Qll), and a fourth on November 28, 2009 (Q12). In 
the following, components Q9, QIO, Qll and Q12 (Fig.[IJ ob- 
served within a distance of 0.5 mas from the core QO are referred 
to as the "jet". 

3.1.1. Kinematics 

We determine the kinematics of individual features from their 
relative positional off'sets with respect to QO. The temporal evo- 



lution of the measured off'sets is plotted in Fig. |3] The compo- 
nent Q8, which was first observed in 2007, is also included for 
the purpose of comparison. 

The motions in the jet of 3C 345 are investigated using the 
R.A., Dec. (x, y) positions of a jet component relative to the 
core component (QO) over the observed periods, fitting them 
separ ately using polynomials of different order (cf . IZensus et al.l 
Il995h and applying the following procedure: 

1 . The first order polynomial is fitted to the trajectory using the 
nonlinear least-squares Marquardt-Levenberg algorithm (as 
implemented by gnuplot). 

2. The second order polynomial is fitted to the trajectory and 
if the improvement of the reduced x~ satisfies the criterion 
described in Section l2!2l the fit is accepted. 

3. The order of the fit polynomial is increased until no improve- 
ment according to the x^ statistics is achieved. 

Using this approach, we find that it is sufficient to represent 
the trajectories of components Q8, Ql 1 and Q12 by linear fits in 
both the X and y directions. In the case of components Q9 and 
QIO, a second order polynomial represents the best fit to the ob- 
served data in the x direction, implying apparent acceleration. 
In the y direction a linear fit is sufficient. The resulting fitted ra- 
dial separations, r(t) - yjx(tY + y(t)-, are drawn in Fig. [3] For 
each component, the fits yield average proper motion and mean 
angular speed {^) and the average direction of motion (<I>). The 
kinematic properties thus derived for the jet components indicate 
that Q9 and QIO underwent a clear phase of apparent accelera- 
tion over a period of 1.5 years and over a distance of ~0.3 mas 
(2pc). No statistically significant acceleration was observed for 
Q8, Ql 1 and Q12. The observed values for (ju) are in the range 
of 0.25 - 0.42 mas yea r'' and for (O) in t he range of -96 to 
-120°, also see Tablel2] |jorstad et alJ(l2005h previouslv reported 
apparent jet component speeds of 0.29 - 0.69 mas year ' and a 
jet position angle of -66 - -95°. Using {jj), the average apparent 
speed yS^pp and the speed in the source frame (deprojected) fi are 
derived through. 



and 



app 



P COS + sin 



(10) 



(11) 



where Z)l is the luminosity distance, z is the redshift and is the 
jet angle to the line of sight (viewing angle). 

As a final step the physical parameters, Doppler factor 6, 
Lorentz factor F and viewing angle of each jet component are 
derived. For the following calculations, we assum e that the jet 
emiss ion is dominated by radiative losses (see e.g.. lJorstad et al.l 
|2005|) and consist of optically thin shocked gas, a - -0.7 (i.e. 
Sv oc R^^-^^^ -v^" , where R is the size of the emitting region in the 
rest-frame of the jet and v the frequency in the observer's frame, 
iMarscher & Geaij|1985h . The variability Doppler factor is then 
derived as. 



^var — 



cAfvai(l +Z)' 



(12) 



where Dl is the luminosity distance and d^ff is the effective an- 
gular size of a spherical region (i.e. the mea sured FXVHM of the 
component multiplied by a factor of 1.8; cf . iPearsonll 1 9991) . The 
flux variability time scale is defined as Afvar = dt/ ln(S 1/82). 
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Fig. 3. Evolution of the radial separations from the core of com- 
ponents Q12, Qll, QIO, Q9 and Q8. Q12-Q9 are related to the 
radio variability since 2008, Q8 is related to a previous radio 
flare observed in 2007. The lines are the results of polynomial 
fits to x(f) and y(f) directions separately. The dashed lines are the 
results of linear fits to x(f) and y{t) to determine ejection epochs 
and speeds for the two cases discussed in Section [3.1.2l 

The value of 5 1 is the measured maximum component flux den- 
sity, while S 2 represents the minimum component flux density 
selected at the time of maximum absolute value of the time 
derivative of the flux density. It is then possible to calculate bulk 
Lorentz factors by combining the derived 5var with jSapp, 



25vi 



(13) 



Combining the derived values for all components, a mean 
Lorentz factor of 12.5 and a mean Doppler factor of 14.4 are 
obtained (see Table|2]i. 

An upper limit for the viewing angle can immediately be 
determined using the average apparent component speed fi^^^ 

and the relation (Smin - This combined with the equa- 

tion 



- arctan 



app 



yS^pp + ^2 _ 1 



(14) 



yields < 5.2°. Using the variability time scale argument 
and the derived values of (5vai from above, the values 0vai for 
each component are obtained, which together h ave a mean 0var 
of_4/7°^g gj, consistent with our upper Umit. iPushkarev et alj 
(I2009h obtained a similar value of = 5.1°, which was deter- 
mined by combining the component speeds at 15 GHz and vari- 
ability Doppler factors derived from single-dish observations at 
37 GHz. Earlier. I Jorstad et al ] (l2005h obtained a smaller viewing 
angle of 2.7°+ 0.9° using VLBI data at 43.2 GHz alone. In the 
following discussion, we adopt = 5°. 

3.1.2. Ejection Epochs 

With the kinematics determined in Section 13.1.11 it is possible 
to estimate the time at which a moving jet feature passes the 
VLBI core at 43 GHz. These passages are referred to as ejection 
epochs. The ejection epoch also marks the time at which a new 
jet feature begins to contribute to the observed radio emission. 



As a first approach to determine the ejection epochs, we as- 
sume that the jet is optically thin all the way to the core and the 
features already travel at the observed average speed while pass- 
ing through the VLBI core. Under these assumptions the ejec- 
tion epoch is determined by back tracing the fitted linear trajec- 
tories. The resulting ejection epochs for all features are listed in 
Table |2] as well. This approach provides reasonable estimates 
of the ejection epochs for non-accelerating features, while for 
the apparently accelerating components Q9 and QIO a different 
approach is required. 

Looking at the radial separations of Q9 and QIO in Fig. [3] the 
acceleration is most evident for the time before 2009.6. An upper 
limit on the ejection epoch can be determined using a linear fit to 
the data points after 2009.5 for Q9 and 2009.4 for QIO (Xupper in 
Table|2]i. Similarly, lower limits, Tio„er, on the ejection epochs of 
Q9 and QIO are obtained by considering only data points before 
2009.5 and 2009.4 respectively. 

3.L3. Jet Intensity Gradient 

In flat spectrum radio quasars (FSRQs), such as 3C 345, the 
component flux decay is c ommonly driven by radiative losses 
dLobanov & ZensusI 119991 : iJorstad et al.ll2005h . which was as- 
sumed in Section [3. 1 . 1 1 without further justification. To test this, 
the maximum component brightness temperature needs to be 
calculated as a measure for the emission intensity of each com- 
ponent, using: 



Tb =1.22- 10' 



' comp 



(l+z) 



dp- ■ v2 



(15) 



where 5 comp is the component flux density in Jansky, z the red- 
shift of the source, (icomp the FWHM size of the circular Gaussian 
in mas and v the observing frequency in GHz. In Fig. |4] the 
brightness temperatures for components Q9, QIO and Qll are 
plotted as a function of the radial separation from the core. 

In the common picture of the shock-in-jet model, a rela- 
tivistic shock propagates down a conical jet, slowly expanding 
adiabatically albeit maintaining shock conditions. In this sce- 
nario the assumption of a power-law electron energy distribution 
{N{E) dE oc £■ ' dE), a power-law magnetic field evolution {B oc 
rr") and a constant jet opening angle with the jet transverse size 
proportional to the distance along the jet (djet 

oc Tjet sin 0) can 

be made. While the shock continues to travel down the jet, it un- 
dergoes three major evolutionary stages d ominated by Compton , 
synchrotron and adiabatic energy losses dMarscher et al.l[T992h . 
From this it follows that the brightness temperature decays as 
a power-law, rb.jet °^ rT^^, where rjet is the distance in the jet at 
which rbjet is measured. The value of e can be derived from 
spectral evolution of radio emission (iLobanov & Zensusll 19991) . 
assuming T\, oc S com^ I {r^^^v^) and the Doppler factor 6 - const. 

For Compton (ec), synchrotron (Es) and adiabatic (fa) losses, 
e are calculated as follows. 



ec = [(i-H5)H-a(i-H l)]/8, 
e, = [A{s + 2) + l>a{s + \)]l6 
fa = [2{2s+l) + 3a{s + \)]l6 



(16) 
(17) 
(18) 



with a typical value of s = 2.0 (corresponding to a synchrotron 
spectral index a - -0.5; S'y oc v^") and a - \ (dominant trans- 
verse magnetic field), - 10/8 = 1.25, - 25/6 ~ 4.17 and 
fa - 19/6 ~ 3.17. For the slopes to be shallower than these de- 
rived values, acceleration resulting in increasing Doppler factors 
of 5 oc r'? may be considered. For a moderate acceleration (i.e. 

jet 
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Table 2. Measured and derived physical parameters for radio emission regions in the inner jet <0.7 mas (4.6 pc). 



Label 


Q8 


Q9 


QIO 


Qll 


Q12 


# 


26 


30 


23 


13 


7 


< yU > (jS^pp) 


0.248 ±0.022 (8.5c) 


0.277 ±0.012 (9.5c) 


0.300 ±0.011 (10c) 


0.427 ±0.020 (15c) 


0.272 ± 0.027 (9.0c) 


/^lower --^ 


- 


0.1260 ±0.0088 


0.165 ±0.016 


- 


- 


upper 


- 


0.458 ±0.089 


0.360 ±0.049 


- 


- 


< <D > 


-100.0° ±2.9° 


-96.3° ±1.1° 


-114.3° ±1.8° 


-120.5° ±3.5° 


-117.8° ±5.4° 


< flower > 


- 


-97.0° ±1.0° 


-109.7° ±3.2° 


- 


- 


< *upper > 


- 


-96.2° ±1.6° 


-112.3° ±2.1° 


- 


- 




2007.62 ±0.12 


2006.91 ±0.11 


2009.055 ±0.025 


2009.544 ±0.018 


2009.651 ±0.042 


"^lower 


- 


2007.834 ±0.083 


2008.674 ±0.073 


- 


- 


"^upper 


- 


2009.00 ±0.16 


2009.11 ±0.11 


- 


- 


Si 


0.73 ±0.12 


3.44 ±0.16 


2.301 ±0.072 


1.489 ±0.086 


2.01 ±0.11 


S2 


0.257 ±0.078 


0.311 ±0.038 


0.816 ±0.044 


0.786 ±0.057 




di 


0.1 14 ±0.038 


0.1115 ±0.0088 


0.1012 ±0.0067 


0.0620 ±0.0094 




d2 


0.244 ± 0.094 


0.232 ±0.036 


0.170±0.016 


0.06 ±0.011 




dt 


1.33 


0.715 


0.385 


0.204 






1.27 ±0.42 


0.297 ±0.016 


0.371 ±0.022 


0.319 ±0.046 




<5var 


5.6 ±2.2 


23.3 ±2.2 


16.9 ± 1.4 


11.9±2.5 




r 

A var 


9.8± 1.4 


13.61 ±0.93 


11.6 ±0.5 


15.1 ± 1.0 




®var 


q 4+2.4 


^•'"-0.23 


9 Q7+0.42 
^•'^'-0.32 






(*) 


no 


yes 


yes 


no 


no 



Notes: # - number of data points; 

< ^ > (jSapp), < yUiower >^ < /"upper > " average compoueut speed in mas year"', in parentheses the apparent speed y6^pp, for Q9 and QIO lower and 
upper limits of the component speeds; 

<<!)>,< Oio„er >> < "tuppcr > - average position angle of the component motion in degree, for Q9 and QIO lower and upper limit; 

Te (Tiower, T„pper) - Ejectiou epoch/timc of zero separation from component labeled QO, determined through linear extrapolation using the average 

component speed and component position at the midpoint, t^^^ = 0.5 • (f^ax - tmm), of the respective dataset, for Q9 and QIO lower and upper 

limits; 

Si, S2 - maximum, minimum (at the time of maximum absolute value of the time derivative) component flux density in units of Jy; 
di, d2 - circular Gaussian modelfit component FWHM in units of mas at the time of Si, S2; 
dt - time between S,^ax and Sn,in in units of years; 
Afvar- - time variability factor in units of years; 

<^vari Tvai-, ©var " derived quantities in order of the variability Doppler and Lorentz factors, and the viewing angle to the line of sight in degrees 
(*) - indicator whether component was accelerating (yes) or not (no) 

The asymptotic standard errors for the fits were determined statistically using an implementation of the nonlinear least-squares 
Marquardt-Levenberg algorithm with the position errors as weights. Errors of the derived parameters are 68% confidence limits determined 
through Monte-Carlo simulations. 



b — 0.1 - 0.2), a value ~ 0.13 is obtained assuming s — 2 and 
with a longitudinal magnetic field (a = 2). 

The brightness temperature gradient along the jet, shown in 
Fig. 3] reveals a possible broken power-law behavior of the jet 
intensity gradient with a break distance of ~0.3 mas and the two 
slopes ei = 0.95 + 0.69 and 62 = 4.11 + 0.85. The value of ei is 
consistent with ec, with the possible indication of a mild change 
in 6, and €2 is consistent with the derived e^. However, at the 
~ Icr level, €2 is also cons istent with e^. Using the sp ectral evo- 
lution of a jet component. iLobanov & ZensusI (Il999h found ev- 
idence for a change from the synchrotron to the adiabatic stage 
at a distance of 1.2-1.5 mas from the core. This suggests that 
the transition from Compton to synchrotron stages is indeed ob- 
served. Future investigations on the spectral evolution of the jet 
are able to provide additional constraints on the transition re- 
gions between the Compton, synchrotron, and adiabatic stages. 

3.2. Radio-y-ray correlation 

The c orrelation be tween radio and optical emission was stud- 
ied in lSchinzel et al . (2011) for identification of y-ray emission 
from 3C345. Here the connection between the synchrotron ra- 
dio emission and y-ray emission is investigated. To quantify a 



1e+13 



1e+12 



1e+11 



Si 1e+10 



1e+09 



1e+08 




Radial Separation (mas) 

Fig. 4. Component brightness temperatures against radial sepa- 
ration from the VLBI core (QO), representing the emission in- 
tensity gradient along the jet. Points with arrows are Icr upper 
limits. Two lines are fitted to the data to determine the power 
law indices e (T^ cc c/rp from the data, the fit from 0. 1 to 0.3 mas 
yields ei = 0.95+0.69 and 0.3 -0.65 mas yields 62 =4.11+0.85. 
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2008 2008.5 2009 2009.5 2010 




2008 2008.5 2009 2009.5 2010 

Epoch (year) 

Fig. 5. top: Fermi LAT 7-day binned y-ray light curve of 3C 345 
for the energy range of 0. l-300GeV. bottom: VLBA 7 mm com- 
ponent flux densities for the modelfitted VLBI core and inner 
jet, represented by up to five circular Gaussian (QO, Q12, Qll, 
QIO, Q9). The component labeled QO is the east-most compo- 
nent (see Fig. [TJ and represents the compact "core" or base of 
the jet. The black curve plots the sum of the flux densities of all 
four components. 



possible connection, the correlation between the y-ray and the 
radio light curve is calculated. For this we used the discr ete cor- 
relation function presented by lEdelson & KrolikI (Il988h . which 
works well for unevenly sampled data. No statistically signifi- 
cant correlation is found, although around 2009.4, a rapid change 
in VLBI flux densities by ~1.5 Jy within only 4 days coincided 
well with a y-ray flare (also see Fig.|5]l. 

The top panel of Fig. [5] plots the weekly averaged y-ray 
light curve of 3C345. From it a seemingly long-term trend is 
seen, most evident between 2008.7 and 2009.4. After 2009.4 
an increase in short-time variability is observed, followed by a 
slightly decreasing trend after 2009.8. The measured y-ray flux 
increases by a factor of 2-3 between 2008.6 and 2009.5. Even 
more evident is a long-term trend in the light curve of the ra- 
dio jet features (bottom panel of Fig. |5]l, following from the sum 
of the flux density of individual jet components. The optical R 
band data presented in Schinzel et al. (201 1) similarly showed a 
rising trend for the same time period. This dataset is added for 
comparison as well. 

To quantify this relation, the light curves are first rescaled 
with the zero flux representing their respective mean value. The 
radio flux density of the jet has a mean value of 3.60 Jy and 
corresponds to the emission from the jet of an apparent size of 
< 2 pc (excluding the core). The core has a mean flux density of 
1.5 Jy and an apparent size of ~ 0.3 pc. Note, the jet flux density 
is by a factor of 2.4 stronger than that of the core. The average 
y-ray flux is 1.8 ■ 10"^ ph cm"^ s"'. Observations at optical R 
band have an average magnitude of 16.6, which corresponds to 
a flux density of 7.7 TO"'* Jy. 

In the final step, the rescaled light curves are de- trended us- 
ing c ubic spline interpolations over 0.4 year bins dPress et alJ 
Il992h . The extracted long-term trends are shown in Fig|6l From 
these it is immediately evident that a similar trend is obtained for 
the radio jet and the y-rays, whereas the core does not show a sig- 
nificant trend. The discrepancy of the first 0.4 year is caused by 
the lack of y-ray monitoring before 2008.6 and a flare immedi- 
ately after 2008.6 (see Fig.|5]l. The trend of the optical emission 



2009 2009.5 2010 




Epoch (year) 

Fig. 6. Long-term trends of the radio jet, radio core flux den- 
sities, optical R band flux (Schinzel et al. 2010) and the y-ray 
flux relative to their respective mean values (radio jet: 3.6 Jy, 
radio core; 1.5 Jy, optical R: 16.6'" (7.7 ■ lO"'* Jy), y-ray; 1.8 ■ 
10"^ ph cm - s '). The trend was extracted fitting cubic splines 
with 0.4 year bins. The relative radio flux density was scaled by 
a factor of 2:1 with respect to the relative y-ray flux, the optical 
flux density was scaled by a factor of 4:1, and the data points 
match the sampling of the original light curves. 



shows a similar behavior compared to the y-ray and radio trends, 
however it is more peaked. Comparing trend amplitudes, the am- 
plitude of the variation in the optical is 4 times greater than in 
the y-rays. The radio shows a factor of 2 higher amplitude of the 
variation than y-rays. Based on the peak values of these trends, 
the radio leads the y-ray trend by 3 1 ^ days and optical leads the 
y-ray trend by l.l^^'^-' days. Altogether this is consistent with an 
almost zero time lag between the observed long-term trends. 

The observed emission from the radio jet showing the match- 
ing trend is most of the time dominated by the jet component 
closest to the core. At the beginning of the monitoring period 
this is Q9 & QIO and by the end of 2010 they are replaced by 
Qll & Q12. This connects the emission from the radio jet with 
that observed at optical and y-ray energies and places the site 
of the underlying multi-wavelength emission within the resolved 
43 GHz radio jet. The very region that was shown to be Compton 
loss dominated in Sec tion l3. 1.31 



4. Discussion 

Monitoring of the fine structure of the jet in 3C 345 reveals a 
wealth of morphological and kinematic features of the flow. In 
particular, compelling evidence for a stationary feature in the 
parsec-scale radio jet of 3C 345 is found. This stationary fea- 
ture is located at a distance of ~0.1 mas (~ 0.7 pc) from the core. 
The kinematics of a large sample of relativistic jets, obtained 
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through the MO J AV E surve^Q, reveale d stationary features as a 
frequent occurence (iLister et al.ll2009l) . In particular, in the jets 
of y-ray blazars, observed at 22 and 43 GHz, stationary features 
within 2 PC (proiec ted) from the core are commonly observed 
jJorstad et al.ll200l"l) . Stationary jet featu res were investigated i n 
more detail in the cases of 3C 390.3 (lArshakian et al.ll20Tol) . 
a nearby double-peaked radio gala xy at z = 0.0561, and the 
Seyfert galaxy 3C 120 (z = 0.033: .Le6n-Tavares et al.l2010l) . For 
these two sources, stationary features were found at distances 
of ~0.4pc (3C 390.3) and ~ 1.3 pc (3C 120) from the core at 
15 GHz. Their physical origin was associated with that of stand- 
ing shocks. A more detailed investigation of the spectral evo- 
lution of a radio flare in CTA 102, a highly polarized quasar at 
z= 1.037, provides supporting evidence for a shock-shock inter- 
action scenario, providing a physical des cription for stationa ry 
features interacting with traveling shocks (iFromm et al.ll201 ll) . 

The moving plasma condensations (components) are most 
likely generated near the base of the jet of 3C 345, but can only 
be reliably monitored after they pass through the location of the 
stationary feature. The nucleus of the jet lies at a distance of 
~7 pc upstream from the 43 GHz VLBI core (Lobanov 1998b; 
ISchinzelll201 ih . After the passage of the stationary feature, the 
superluminal components in 3C 345 undergoe an apparent ac- 
celeration from ~5 to ~ 15 c over a distance of 0.3 mas in ex- 
tent (~23pc deprojected). This acceleration seems to be a com- 
monly observed feature of many objects in the MOJAVE sam- 
ple, where a statistically sign ificant tendency for acceleration 
in the base of jets is fou nd (iHoman et al.l l2009h . Previously, 
ILobanov & Rolandl ( l2005h discussed the observed acceleration 
of a different component in 3C 345, testing the case of a sub- 
stantially curved three-dimensional path, causing the observed 
dramatic change of the apparent component speed. They con- 
cluded that the observed acceleration is not likely to be a geomet- 
rical eff'ect, but should reflect physical acceleration of plasma. 
The possibility of intrinsic acceleration of jet feature s on par- 
sec scales was discussed bv lVlahakis & Konigll (l2004 . There it 
was specifically argued that the acceleration observed in 3C 345 
is not purely hydrodynamic, but can be attributed to magnetic 
driving. They showed that the observed acceleration is consis- 
tent with an acceleration from F ~ 5toF ~ 10 over a linear 
distance of ~3 to ~20pc. The results presented here are consis- 
tent with these findings. 

The times of passage of knots QIO, Qll, and Q12 through 
the 43 GHz VLBI core coincide with short-term increases in the 
■y-ray flux (see Fig.|5] with times listed in Table |2]i. In a similar 
case, PKS 1510-089 showed, in addition to the appearance of 
a new VLBI component, a rapid rotation in the optical linear 
polarization just prior the components passage through the core. 
This was explained by the feature following a spiral path through 
a tor oidal field until it cro sses a standing shock in the 43 GHz 
core dMarscher et al.ll2010l) . 

A y-ray flare centered on 2009.4 may be related to the bright- 
ening of radio structure at distances of 0.12 - 0.20 mas from the 
core (0.8 - 1.3 pc), corresponding to a deprojected distance of 
~9-15pc. Interestingly, this includes the region for which the 
presence of a standing shock is suggested by the observations 
presented in Sec. 13.1.11 Radio flares in the comp act jet down- 
stream of the VLBI core are known to occur (e .g., IWehrle et al.l 
1200 ll) . Furthermore. Ijorstad et all (1200 ll 12010) have found that 
y-ray flares often occur after superluminal knots have separated 
from the VLBI core. The results presented here support the 
inference that y-ray emission still occurs parsecs downstream 
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of the central engine, contrary to conclusions based on theo- 
retical arguments (e.g.. lTavecchio et"ani2010HPoutanen & StemI 
2010). These arguments place the y-ray emission zones within 
the broad-line region (BLR) at distances of < 1 pc from the cen- 
tral engine. Within this region a high numb er of optical-UV pho- 
tons are produced (e.g. lKaspi et al.ll2007l) . providing ideal con- 
ditions for inverse Compton scattering with relativistic electr ons 
from the jet (e.g. [Arshakian et al. 201 1; Dermer et an il997l) . In 
addition, the observed variability time scales place constraints 
on the size of the emitting region due to causality arguments. 
Usually this size is related to the jet cross section, implying that 
the emission region is located close to the black hole. However, 
these standard assumptions do not necessarily need to be true. 

In Section l3T2l we have shown that long-term trends observed 
at y-ray and optical energies correlate with those of the radio 
emission at VLBI scales, which imply a common emission re- 
gion. This emission region is related to the inner jet, extending 
over a distance of up to 23 pc from the VLBI core, which is at 
a distance of »1 pc from the central engine, well beyond the 
BLR. The mechanism for the emission at radio and optical wave- 
lengths is therefore dominated by the synchrotron emission from 
the relativistic jet itself. The simultaneity of y-ray emission im- 
plies inverse Compton upscattering of synchrotron photons pro- 
duced in the jet. The relative increase of optical and radio emis- 
sion compared to y-ray emission is a factor of 2 higher in the 
radio and even a factor of 4 higher in the optical emission, which 
is consistent with the mechanism of synchrotron self-Compton. 
An external contribution of "se ed" photons, as propo sed by ex- 
ternal Compton scenarios (e.g. iTavecchio et aljl2010l) is not re- 
quired to produce the amount of the observed y-ray emission. 
Additionally, we have found that the properties of the emis- 
sion in this region are consistent with the Compton-loss stage 
of a shock (Section r3.1.3b . Altogether, this contradicts the com- 
mon theoretical high energy production scenario described in the 
previous par agraph. Recent obse rvations of TeV emission from 
FSRQs (e.g. lAleksic et alJl201 ll) should not have been possible 
if the high energy emission site were within the BLR, which is 
expected to be opaque to y-rays at TeV energies due to yy inter- 
actions. 

The observed kinematics of the jet is now combined with the 
rise times of the radio, optical, and y-ray emission. The com- 
ponent Q9 is the first new component identified in 2008 and is 
most likely the emission region responsible for the initial flux 
increase. Using the long-term trends extracted, the distance at 
which Q9 was located during the onset of the flaring period in 
2008 is ~0.05 mas from the VLBI core. This marks the region 
between the stationary feature and the core at 43 GHz. Instead 
of using the combined flux trend that is related to the sum of 
the radio emission regions, we calculate the time between the 
peak brightness of Q9 (2009.40) and the first y-ray flare ob- 
served (2008.73). If we use the average component speed and 
assume no acceleration, then the closest Q9 could have been dur- 
ing the time of the first observed y-ray flare is consistent with the 
position of the radio core at 43 GHz. This finding strongly sup- 
ports the interpretation of the parsec-scale jet being responsible 
for driving the o b served high energy emission. 

lAgudo et alj (1201 Ih have reported on the location of a y- 
ray flare in OJ287, which was placed at a distance of more 
than 14pc from the central engine. They proposed a model for 
the multi-wavelength emission where the high energy emission 
could be explained by the synchrotron self-Compton process or 
inverse Compton scattering of infrared radiation from a hot dusty 
torus, although it was concluded that the hot dusty torus scenario 
is less likely. In the case of 3C 345 it seems to be clear that mul- 
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tiple compact emission features are responsible for the observed 
variability and what seems to be a reasonable model for the case 
of OJ 287, could apply for 3C 345 as well. A new plasma distur- 
bance passes through a first conical shock in the core and pro- 
duces a fast rise in high energy emission; as it continues to prop- 
agate down the jet it continues to produce high energy emission. 

5. Summary & Conclusions 

Based on 32 VLBA observations of 3C 345 at 43.2 GHz, we have 
investigated the structure and evolution of the radio emission 
and related this to variations of the y-ray and optical emission 
observed 2008-2010 by Ferm;-LAT and a number of optical ob- 
servatories. 

We have identified and analyzed four new moving emission 
features (jet components) in the radio jet of 3C345. These re- 
gions are found to move at apparent speeds of 9-15 c, with the 
corresponding Doppler and Lorentz factors derived to be in the 
ranges of 12-23 and 12-15, respectively. The kinematic data 
strongly favor a viewing angle of 4.7° between the jet axis and 
the line of sight. 

We have presented evidence for the y-ray emission to be pro- 
duced not in a compact region near the central engine of the 
AGN, but in the Compton-loss dominated zone of the parsec- 
scale jet extending up to r 0.3 mas, corresponding to a depro- 
jected linear extent of x 23 pc (accounting for the source distance 
and jet orientation). This zone is further marked by a break in the 
evolution of the brightness temperature, Tb, of the radio emis- 
sion, with Tb oc ^-(0-95 ±0.69) ^ < 0.3 mas and Tb oc ^-(4.ii±o.85) 

at larger separations. Ejections of new superluminally moving 
and apparently accelerating features in the jet are linked to the 
flaring component of the y-ray and optical emission. 

These findings favor the synchrotron self-Compton mecha- 
nism of the high-energy emission production, while questioning 
the entire class of models that place the high energy emission 
site within 1 pc from the central engine of AGN, with external 
seed photons for inverse Compton scattering from the accretion 
disk or broad-line-region. They also imply that a significant part 
of the y-ray emission is generated by highly-relativistic elec- 
trons propagating at a large bulk speed in the jet. In the context 
of newly emerging results, stimulated by Fermi data, more de- 
tailed analytical and numerical descriptions of such a scenario 
are clearly needed in order to explain the observed connection 
between radio and y-ray variability. At the same time, contin- 
ued monitoring and more densely sampled VLBI observations 
of well studied, bright blazar radio jets represent an essential re- 
quirement for improving the degree of detail and statistical accu- 
racy of the correlation reported and further improving the spatial 
locahzation of individual flares in relativistic jets. 

Acknowledgements. We thank Nicola Mai'chili for providing us with his imple- 
mentation of the cubic spline de-trending algorithm. We thank the anonymous 
referee for valuable comments. Frank Schinzel was supported for this research 
through a stipend from the International Max-Planck Research School (IMPRS) 
for Astronomy and Astrophysics at the Universities of Bonn and Cologne. We 
thank NASA for support under FERMI grant GSFC #21078/FERMI08-005I. 
The research at Boston University was supported in part by NASA through Fermi 
Guest Investigator grants NNX08AJ64G, NNX08AU02G, NNX08AV61G, and 
NNX08AV65G, and by National Science Foundation grant AST-0907893. The 
National Radio Astronomy Observatory is a facility of the National Science 
Foundation operated under cooperative agreement by Associated Universities, 
Inc. This research has made use of the NASA/IPAC Extragalactic Database 
(NED) which is operated by the Jet Propulsion Laboratory, California Institute 
of Technology, under contract with the National Aeronautics and Space 
Administration. This research has made use of NASA's Astrophysics Data 
System. 



References 

Abdo, A. A., Ackermann, M., Ajello, M., et al. 2010, ApJS, 188, 405 

Abdo, A. A., Ackermann, M., Ajello, M., et al. 2009, ApJ, 700, 597 

Agudo, I., Jorstad, S. G., Marscher, A. P, et al. 2011, ApJ, 726, L13 

Aleksic, J., Antonelli, L. A., Antoranz, P, et al. 2011, ApJ, 730, L8 

Aller, H. D., AUer, M. E, & Hughes, P. A. 1996, in Astronomical Society of 

the Pacific Conference Series, Vol. 110, Blazai' Continuum Variability, ed. 

H. R. Miller, J. R. Webb, & J. C. Noble, 208 
Arshakian, T. G., Leon-Tavai'es, J., Boettcher, M., et al. 2011, eprint 

arXiv:l 104.4946 

Arshakian, T. G., Leon-Tavares, J., Lobanov, A. R, et al. 2010, MNRAS, 401, 
1231 

Atwood, W. B., Abdo, A. A., Ackermann, M., et al. 2009, ApJ, 697, 1071 
Avni, Y. 1976, ApJ, 210, 642 

Baath, L. B., Rogers, A. E. E., Inoue, M., et al. 1992, A&A, 257, 31 
Babadzhanyants, M. K. & Belokon, E. T. 1984, Astrophysics, 20, 461 
Biretta, J. A., Moore, R. L., & Cohen, M. H. 1986, ApJ, 308, 93 
Blazejowski, M., Siemiginowska, A., Sikora, M., Moderski, R., & Bechtold, J. 

2004, ApJ, 600, L27 
Bottcher, M., Reimer, A., & Marscher, A. P 2009, ApJ, 703, 1168 
Casandjian, J.-M. & Grenier, I. A. 2008, A&A, 489, 849 
Dermer, C. D., Finke, J. D., Krug, H., & Bottcher, M. 2009, ApJ, 692, 32 
Dermer, C. D., Stumer, S. J., & Schlickeiser, R. 1997, ApJS, 109, 103 
Edelson, R. A. & Krohk, J. H. 1988, ApJ, 333, 646 

Fomalont, E. B. 1999, in Astronomical Society of the Pacific Conference Series, 
Vol. 180, Synthesis Imaging in Radio Astronomy II, ed. G. B. Taylor, 
C. L. CarilH, & R. A. Perley, 301 

Fromm, C. M., Perucho, M., Ros, E., et al. 2011, A&A, 531, A95+ 

Hirotani, K. 2005, ApJ, 619, 73 

Homan, D. C, Kadler, M., Kellermann, K. I., et al. 2009, ApJ, 706, 1253 
Jorstad, S. G., Mai'scher, A. P, Lister, M. L., et al. 2005, AJ, 130, 1418 
Jorstad, S. G., Mai'scher, A. P, Mattox, J. R., et al. 2001, ApJS, 134, 181 
Jorstad, S. G., Marscher, A. P, Smith, P A., Larionov, V. M., & I., A. 2010, in 

Fermi meets Jandky - AGN in Radio and Gamma Rays, ed. T. Savolainen 
Jorstad, S. G., Mai'scher, A. P, Stevens, J. A., et al. 2007, AJ, 134, 799 
Kaspi, S., Brandt, W. N., Maoz, D., et al. 2007, ApJ, 659, 997 
Kidger, M. & Takalo, L. 1990, A&A, 239, L9 

Klare, J. 2004, PhD thesis, Rheinische Fiiedrich-Wilhelms-Universitat Bonn, 

Bonn, Germany, (2004) 
Larionov, V. M., Villata, M., Raiteri, C. M., et al. 2009, The Astronomer's 

Telegram, 2222, 1 

Leon-Tavares, J., Lobanov, A. P, Chavushyan, V. H., et al. 2010, ApJ, 715, 355 
Lister, M. L., Cohen, M. H., Homan, D. C, et al. 2009, AJ, 138, 1874 
Lobanov, A. P. 1996, PhD thesis, , New Mexico Institute of Mining & 

Technology, Socorro, NM, US, (1996) 
Lobanov, A. P 1998a, A&AS, 132, 261 
Lobanov, A. P 1998b, A&A, 330, 79 
Lobanov, A. P. 2005, eprint arXiv:astro-ph/0503225 
Lobanov, A. P & Roland, J. 2005, A&A, 431, 831 
Lobanov, A. P & Zensus, J. A. 1999, ApJ, 521, 509 

Longair, S. M. 201 1, High Energy Astrophysics (Cambridge University Press) 
Marscher, A. P & Gear, W. K. 1985, ApJ, 298, 114 

Marscher, A. P, Gear, W. K., & Travis, J. P. 1992, in Variability of Blazars, ed. 

E. Valtaoja & M. Valtonen, 85 
Marscher, A. P, Jorstad, S. G., Larionov, V. M., et al. 2010, ApJ, 710, L126 
Marziani, P., Sulentic, J. W., Dultzin-Hacyan, D., Calvani, M., & Moles, M. 

1996, ApJS, 104, 37 
Pearson, T. J. 1999, in Astronomical Society of the Pacific Conference Series, 

Vol. 180, Synthesis Imaging in Radio Astronomy II, ed. G. B. Taylor, 

C. L. CarilH, & R. A. Perley, 335 
Poutanen, J. & Stern, B. 2010, ApJ, 717, LI 18 

Press, W. H., Teukolsky, S. A., Vetteriing, W. T., & Flannery, B. P 1992, 
Numerical recipes in FORTRAN. The art of scientific computing, ed. Press, 
W. H., Teukolsky, S. A., Vetteriing, W. T, & Flannery, B. P 

Pushkarev, A. B., Kovalev, Y. Y, Lister, M. L., & Savolainen, T. 2009, A&A, 
507, L33 

Ros, E., Zensus, J. A., & Lobanov, A. P 2000, A&A, 354, 55 
Schinzel, F. K. 201 1, PhD thesis. University of Cologne 

Schinzel, F. K., et al., & for the Fermi-LAT collaboration. 201 1, A&A, accepted, 

to appear in preprints and on astro-ph 
Shepherd, M. C, Pearson, T. J., & Taylor, G. B. 1995, in Bulletin of 

the American Astronomical Society, Vol. 27, Bulletin of the American 

Astronomical Society, ed. B. J. Butler & D. O. Muhleman, 903 
Tavecchio, E, Ghisellini, G., Bonnoli, G., & Ghirlanda, G. 2010, MNRAS, 405, 

L94 

Teraesranta, H., Tornikoski, M., Mujunen, A., et al. 1998, A&AS, 132, 305 
Thompson, D. J., Bertsch, D. L., Eichtel, C. E., et al. 1993, ApJS, 86, 629 



10 



Schinzel et al.: Relativistic Outflow Drives y-Ray Emission in 3C 345 



Unwin, S. C, Wehrle, A. E., Lobanov. A. P., et al. 1997, ApJ, 480, 596 

Valtaoja, E. & Teraesranta, H. 1996, A&AS, 120, C491 

Vlahakis, N. & Konigl, A. 2004, ApJ, 605, 656 

Wehrle, A. E., Finer, B. G., Unwin, S. C, et al. 2001, ApJS, 133, 297 

Zensus, J. A., Cohen, M. H., & Unwin, S. C. 1995, ApJ, 443, 35 



